Interference of Bose-Einstein condensates in momentum space 
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We suggest an experiment to investigate the linear su- 
perposition of two spatially separated Bose-Einstein conden- 
sates. Due to the coherent combination of the two wave func- 
tions, the dynamic structure factor, measurable through in- 
elastic photon scattering at high momentum transfer q, is 
predicted to exhibit interference fringes with frequency period 
Av = q/md where d is the distance between the condensates. 
We show that the coherent configuration corresponds to an 
eigenstate of the physical observable measured in the experi- 
ment and that the relative phase of the condensates is hence 
created through the measurement process. 

PACS numbers: 03.65.-w, 05.30.J, 32.80 -t, 67.40.-w 



Dilute Bose-Einstein condensed gases behave like clas- 
sical matter waves. This striking feature has been di- 
rectly confirmed by several recent experiments In 
particular, very clean interference patterns generated by 
two overlapping condensates have been observed through 
absorption imaging techniques Q . Interference phenom- 
ena produced by matter waves are key features underly- 
ing the quantum mechanical behaviour of matter, so it 
is of considerable interest to understand what is the new 
role played by Bose-Einstein condensation. 

Bose-Einstein condensed gases can be regarded as clas- 
sical objects because, according to Bogoliubov prescrip- 
tion, the corresponding field operator can be replaced by 
a classical field, resembling the classical limit of quan- 
tum electrodynamics. Differently from the case of the 
electromagnetic field which is governed by the Maxwell 
equations, the field associated with a Bose-Einstein con- 
densate obeys equations of quantum nature which re- 
duce, for dilute and cold gases, to the Gross-Pitaevskii 
equation 0] 
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where V ex t is the external potential confining the gas, and 
g = 4Trh 2 a/m is the interaction coupling constant, fixed 
by the s-wave scattering length a. The field 'J has the 
meaning of an order parameter and is often called the 
"wave function of the condensate". As a consequence of 
the quantum nature of equation (Q), key features of the 
field 'J', like for example interference patterns, depend 
explicitly on the value of the Planck constant. 

The fact that two overlapping condensates behave as 
coherent matter waves giving rise to interference is not 



however obvious. In a similar context Anderson || raised 
the intriguing question "do two superfluids which have 
never "seen" one another possess a definitive phase?". 
This question has been the object of theoretical specula- 
tions 0| and has been more recently reconsidered |^ [l^] 
after the experimental realization of BEC in trapped 
atomic gases |13| . The point of view shared by most au- 
thors is that the relative phase between two condensates 
is "created" during the measurement. In other words, 
even if the initial configuration is not coherent and the 
relative phase between the condensates is not fixed, one 
can still observe interference in a single realization of the 
experiment. This opens new interesting perspectives in 
the field of quantum measurement in macroscopic sys- 
tems. However, measuring fringe patterns in the density 
profiles requires overlapping of the condensates in coor- 
dinate space and one cannot exclude the possibility that 
interactions among atoms block the relative phase before 
measurement [ ^0|JT^ ]. 

In this paper we propose an alternative way to investi- 
gate interference and coherence effects, by exploring the 
behaviour of the condensate in momentum rather than 
in coordinate space. Our proposal is stimulated by the 
recent experiment of p"5[ where, by measuring the dy- 
namic structure factor at high momentum transfer via 
stimulated two-photon Bragg scattering, it was possible 
to prove that a single condensate does not exhibit phase 
fluctuations, i.e "it does not consist of smaller quasi- 
condensates with random relative phase" p!q |. In fig. 1 
we show a useful configuration where two parallel trapped 
condensates are located at distance d along the x axis and 
are described by the order parameters ^ a and 'Iv The 
condensates have no overlap in coordinate space but can 
exhibit interference in momentum space. Geometries of 
this type are now becoming available via optical confine- 
ment techniques. In the presence of coherence the or- 
der parameter of the whole system is given by the linear 
combination ty c = ^ a + e l ^^>i,- The corresponding many- 
body function is (c^) N |> where c 1 is the particle creation 
operator in the state \I/ C and | > is the vacuum of parti- 
cles. In the following we assume, for simplicity, that the 
two condensates have the same shape, contain the same 
average number of atoms (N a = Nf, = N/2) and that the 
potential separating the two condensates is large enough 
to exclude overlap between the two wave functions. Un- 
der these conditions we can write ^(r) = ^o( r + d/2) 
and 'l'b(r) = ^o(r — d/2) where is the order param- 
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eter of a single condensate obeying the Gross-Pitaevskii 
equation ([j]) and normalized to J dr | Wo | 2 = N/2. In 
momentum space the order parameter takes the form 

*(p) = e- ip * d/2h ^ {p) + e^ +fed / 2ft )* (p) (2) 

where * (p) = (2?:hy 3/2 J dre(~ip ■ r/7L)* (r). . 

Let us calculate the average value of the momen- 
tum distribution operator n(p) = ^^(p)*(p), where 
'J'(p) is the field operator in momentum representation. 
For the coherent configuration (j^) one finds the result 
n(p) =< n(p) >= 2(1 + cos(p x d/K + 0))no(p) which ex- 
hibits interference fringes with period Ap x — 2nh/d [ [l6|| . 
In this equation Uo(p) =| v £'o(p) | 2 is the momentum 
distribution of each condensate. 

In the absence of coherence the many-body wave func- 
tion has instead the form (<v) Na Qr) b |> where and 

are the particle creation operators in the states and 
Wb respectively. This state corresponds to two indepen- 
dent condensates with fixed number of atoms iV a and Nb 
respectively, and the average of the momentum distribu- 
tion operator takes the value n(p) = 2no(p) which, as 
expected, does not exhibit interference. 

The momentum distribution of the condensate can be 
investigated experimentally by measuring the dynamic 
structure factor at high energy and momentum transfer. 
A first useful description is provided by impulse approx- 
imation |l7| ] 



m f 

S(q,E) = — J n(Y,p y ,p z )dpydp z 
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which relates the dynamic structure factor to the so 
called longitudinal momentum distribution v{p x ) — 
J n(p x ,Py,Pz)dpydp z . In @ E and q are the energy and 
momentum transferred by the photon to the system and 

2 

Y = ^(E — t^-) is the relevant scaling variable of the 
problem fjgft . The vector q has been taken along the x 
direction. A remarkable feature of impulse approxima- 
tion is that the quantity qS(q, E) depends on q and E 
only through the scaling variable Y. Impulse approxi- 
mation assumes that the system, after scattering with 
the photon, can be described in terms of a scattered 
atom propagating with momentum p + q and (N — 1) 
atoms remaining in the unperturbed configuration. This 
approximation has been extensively used to analyze the 
momentum distribution of varius classical and quantum 
systems, including liquids and solids through deep inelas- 
tic neutron scattering. In particular it has been employed 
to extract the condensate fraction of superfluid 4 He [[l9| . 
The impulse approximation is accurate if one can ignore 
final state interaction effects which are responsible for 
both a shift of the peak energy with respect to the free 
recoil value E r — q 2 /2m and for a broadening of the 
function S(q, E). First measurements and theoretical es- 
timates of these effects in a cold trapped Bose gas have 



been presented in |Ts|| . For momentum transfers signif- 
icantly smaller than the inverse of the scattering length 
one can safely use Bogoliubov theory. The relative shift 
of the peak is a small effect, fixed by the ratio ii/E r where 
ii = gn(0) is the chemical potential and n(0) is the cen- 
tral density of the gas. Typical values in the experiment 
of (Tj| are E r ~ 20 — 100/x, depending on the density of 
the sample. The broadening due to interactions, of the 
order of the chemical potential, should be compared with 
the Doppler broadening AE^oppier ~ hq/mR x , fixed by 
the width h/R x of the momentum ditribution of the con- 
densate and increasing linearly with q. Here R x is the 
radius of the condensate in the x direction, determined 
by the Thomas-Fermi relation ii = muj 2 R 2 /2, where us x 
is the radial frequency of the harmonic potential trap- 
ping each condensate. The condition derived in [fL5| for 
the Doppler effect being the leading effect can be put in 
the form E r /[i > O.Q2{ii/hio x ) 2 . Since the ratio ii/fous x is 
large for high density samples, this condition can become 
rather severe. For example, in the experiment of |jl5[ it 
is well satisfied only for the low density samples. 

A striking feature predicted by (||) is that, in the pres- 
ence of coherence, also the dynamic structure factor ex- 
hibits interference. In fact, inserting the corresponding 
result for the momentum distribution into ml), one finds 



S(q, E) = 2[1 + cos(Yd/h + (f>)]S {q, E), 
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where So(q, E) is the dynamic form factor relative to each 
condensate. The fringes have a period 
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and their position is fixed by the value of the relative 
phase (f> which can be consequently measured. Notice 
that the ratio between the Doppler width of S(q, E) and 
the distance between two fringes scales as ~ d/R x . A 
smoothing of the interference signal is expected to oc- 
cur if one includes final state interaction effects. It is 
worth noticing that result (J4j) holds only if the energy 
of the scattered atoms, of the order of the recoil energy 
~ q 2 /2m, is larger than the height of the barrier sepa- 
rating the two condensates |^0| . In fact only in this case 
will the reflection of the scattered atom from the barrier 
be negligible and the dynamic structure factor S(q, E) be 
sensitive to the phase <f> of the two condensates. Actually 
the scattered atoms provide the physical coupling be- 
tween the two condensates. In the absence of coherence, 
or if the recoil energy is smaller than the height of the 
barrier, one has instead S(q,E) = 2So(q,E). In fig. 2 we 
show a typical prediction for the dynamic structure func- 
tion in the absence (dashed line) and in the presence of 
coherence (full line) between the two condensates. The 
curves have been calculated for a gas of sodium atoms 
trapped in two identical cigar shaped harmonic traps 
with frequencies v z = 10Hz, v±_ = 100Hz and central 
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density equal to 0.5 x 10 14 cto~ 3 . The distance between 
the centers of the trap is d — AR X ~ 33/im and the mo- 
mentum transfer q is equal to 21.3fi(/xm) _1 . With such 
parameters the broadening of the dynamic structure fac- 
tor due to final state interactions can be shown to be 
smaller than the period of interference fringes. The visi- 
bility of fringes would become even better by increasing 
the momentum transfer q and/or reducing the density of 
the sample. 

It is interesting to compare the interference fringes 
exhibited by the momentum distribution with the ones 
characterizing the density of two expanding and overlap- 
ping condensates jD. A simple estimate of the density 
fringes is obtained by neglecting interaction effects be- 
tween the two condensates during the expansion. This is 
a good approximation if the two condensates are initially 
well separated in space. Let us write the order parameter 
of each expanding condensate in the form Wo — \/ r M)C lx 
where no is the density and \ is the phase. The in- 
terference patterns of the density n(r) =| ^(r) | 2 are 
fixed by the phase difference x( r + d/2) — x(r — d/2). 
In the geometry of fig. 1 this difference approaches very 
rapidly the asymptotic value mxd/ht corresponding to 
a relative velocity v — d/t of the two condensates in 
the x direction ||l],|2^]. In the presence of coherence 
the density profile of the overlapping clouds takes the 
form n(r, t) = n a + rij + 1^Jn a n\, cos(xdm/ht + (f>) where 
n a = n o{ r + d/2, t) and rib = no(r — d/2, t). The result- 
ing fringes are straight lines orthogonal to the a;-axis and 
have wave length A = 2itfit/md. Such fringes have been 
directly observed in the experiment of [ . 

With respect to the experiment of II], based on mea- 
surements of the density of two overlapping condensates, 
the measurement of the dynamic structure factor at high 
momentum transfer should exhibit interference even if 
the condensates are separated in space. This raises the 
question of what happens to S(q, E) if the two conden- 
sates have never "seen" one another, i.e. if they are not in 
a relatively coherent state before measurement. Accord- 
ing to quantum measurement theory, the system, after 
measurement, jumps into an eigenstate of the measured 
observable. So it is important to understand what are the 
eigenvalues of the momentum density operator which is 
the physical observable measured in a single realization 
of the experiment. This should be distinguished from the 
average value taken on several realizations of the experi- 
ment. 

According to the Bogoliubov prescription, a fully Bose- 
Einstein condensed state is not only eigenstate of the field 
operator W, but also of the density as well as of the mo- 
mentum density operators. Of course this prescription 
does not apply to configurations exhibiting fragmenta- 
tion of Bose-Einstein condensation, as happens in the 
case of two independent condensates. Furthermore it 
ignores microscopic fluctuations arising from the non- 
commutativity of W and $K For the above reasons it 



is interesting to discuss in a more quantitative way what 
are the conditions of applicability of the Bogoliubov pre- 
scription. For a proper discussion of the problem it is 
crucial to consider macroscopic coarse grained averages 
of the physical observables. This averaging takes into ac- 
count the finite resolution of the experimental apparatus. 
On the other hand integrating the signal can be crucial in 
order to produce visible interference patterns. From the 
theoretical side one can show that only by taking these 
averages will the fluctuations of the physical observables 
become negligible. Let us discuss in details the problem 
of the momentum distribution which is the main object 
of the present work. A similar discussion can be repeated 
for the fluctuations of the density operator in the context 
of experiments where two expanding condensates overlap 
in coordinate space and are then imaged [Q. Since the 
momentum distribution enters the relevant expression (|^) 
integrated with respect to p y and p z , one only needs to 
consider the coarse grained average along the x direction 
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where n(p) = ^(p)^^) i s the momentum distribution 
operator and, for simplicity, we have chosen a gaussian 
convolution. The fluctuations of the operator Pp(px) are 
easily calculated for a coherent configuration. After or- 
dering the field operators one finds the result 
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where the term in the right hand side arises from the 
non commutativity of the field operators W and WL This 
equation explicitly shows that the fluctuations diverge 
when [3 — > 0. Viceversa they become negligible in the 
macroscopic limit f3v(p x ) 3> 1, i.e. if the number of 
atoms with momentum between p x and p x + (3 is suffi- 
ciently large. This condition is easily satisfied for actual 
Bose-Einstein condensates and is compatible with the re- 
quest that P be smaller than the distance Ap x = 2Trh/d 
between two consecutive fringes of the momentum distri- 
bution. As a consequence of the strong quenching of the 
fluctuations, the coherent state can be considered, with 
proper accuracy, an eigenstate of the operator vp{p x ). 
The situation is very different if one instead considers 
two independent condensates occupying, respectively, the 
single particle wave functions W a and ^b- In fact in this 
case the randomness of the relative phase produces the 
additional contribution < Vp(p x ) > 2 /2 to the fluctua- 
tions of P/3(p x ) which become macroscopically large. 

The above discussion reveals that a state made of 
two independent condensates is not an eigenstate of the 
macroscopic observable (^). As a consequence of mea- 
surement, the system will jump into a coherent configura- 
tion and will exhibit interference patterns. Furthermore, 
if the experiment is non destructive the two condensates 
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will remain trapped and separated in space also after the 
measurement. The possibility of observing interference in 
a single realization of the experiment is a peculiar con- 
sequence of Bose-Einstein condensation. Of course one 
should carry out the measurement within times shorter 
than the phase decoherence times Different real- 

izations of the experiment on independent condensates 
would give rise to different values of the phase and con- 
sequently to strong fluctuations in the measured signal. 

In conclusion we have pointed out the occurrence of in- 
terference phenomena in momentum space exhibited by 
Bose-Einstein condensates separated in coordinate space. 
The interference patterns should be visible in the dy- 
namic structure factor measured in photon scattering 
experiments at high momentum transfers. Such exper- 
iments would open new intriguing perspectives in the 
study of quantum superposition and non locality phe- 
nomena in macroscopic systems. 
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FIG. 1. Typical configuration of two sepa- 
rated Bose-Einstein condensates, confined by parallel cigar 
shaped traps. 



4 




-4-2 2 4 

v— v r [kHz] 

FIG. 2. Dynamic structure factor calculated with (full line) 
and without (dashed line) coherence between the two conden- 
sates. The dynamic structure factor was calculated in impulse 
approximation, using the Thomas-Fermi limit for the ground 
state momentum distribution (see text). The relative phase 
was taken equal to zero, and v = E/2nh. 
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